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Abstract 



(— I , A new functional for simplicial surfaces is suggested. It is invariant 

with respect to Mobius transformations and is a discrete analogue of 
the Willmore functional. Minima of this functional are investigated. 
As an application a bending energy for discrete thin-shells is derived. 
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1 Introduction 

In the variational description of surfaces the following functionals are of 
primary importance: 

• The area A = J dA, where dA is the area element, is preserved by 
isometries. 

• The total Gaussian curvature Q = J KdA, where K is the Gaussian 
curvature, is topological invariant. 



• 



The total mean curvature M = J HdA, where H is the mean curva- 
ture, depends on external geometry of the surface. 

• The Willmore energy W = / H^dA is invariant with respect to Mobius 
transformations. 

Geometric discretizations of the first three functionals for simplicial sur- 
faces are well known. For the area functional it is obvious. The local Gaus- 
sian curvature at a vertex v is defined as the angle defect 

G{v) = 27:-Y,ai, 

i 

where ckj are the angles of all triangles (see Fig. EJ at vertex v. The total 
Gaussian curvature is a sum over all vertices G = ^^ G{v). The local mean 
curvature at an edge e is defined as 

M(e) = le, 

where / is the length of the edge and 6 is the angle between the normals 
to the adjacent faces at e (see Fig. ^. The total mean curvature is the 
sum over all edges M = 'Ylie^i^)- These discrete functionals possess the 
geometric symmetries of the smooth functionals mentioned above. 

Until recently a geometric discretization of the Willmore functional was 
missing. In this paper we introduce a Mobius invariant energy for simplicial 
surfaces and show that it should be treated as a discrete Willmore energy. 

2 Conformal Energy 

Let S be a simplicial surface in 3-dimensional Euclidean space with the set 
of vertices V , edges E and (triangular) faces F. We define a conformal 
energy for simplicial surfaces using circumcircles of their faces. Each (inter- 
nal) edge e G i? is incident to two triangles. Consistent orientation of the 
triangles naturally induces an orientation of the corresponding circumcir- 
cles. Let /3(e) be the external intersection angle of the circumcircles of the 
triangles sharing e, which is the angle between the tangent vectors of the 



oriented circumcircles. Now we are ready to define a new energy functional 
for simplicial surfaces. 

Definition 1 Local conformal (discrete Willmore) energy at a vertex v is 
given by the sum over all incident edges 



W{v) = ^/3(e)-27r. 
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The conformal (discrete Willmore) energy of the simplicial surface is the 
sum over all vertices 



vev eeE 

where \V\is the number of vertices of S. 





Figure 1: Definition of conformal (discrete Willmore) energy 

Fig^presents two neighboring circles with their external intersection an- 
gle Pi as well as a view "from the top" at a vertex v showing all n circumcir- 
cles passing through v with the corresponding intersection angles /?i, . . . , /3„. 
For simplicity we will consider only simplicial surfaces without boundary. 

Note that the energy we defined is obviously invariant with respect to 
Mobius transformations. This invariance is an important property of the 
classical Willmore energy defined for smooth surfaces (see below). 

Note also that the changing of the orientation of both circles preserves 
the angle /3(e). This shows that the energy W{S) is well defined for non- 
oriented simplicial surfaces as well. 

The star S{v) of the vertex v is the subcomplex of S comprised by the 
triangles incident with v. The vertices of S{v) are v and all its neighbors. 
We call S{v) convex if for any its face / G F{S{v)) the star S{v) lies to one 
side of the plane of F and strictly convex if the intersection of S{v) with the 
plane of / is / itself. 



Proposition 2 The conformal energy is non-negative 

W{v) > 0, 

and vanishes if and only if the star S{v) is convex and all its vertices lie on 
a common sphere. 

The proof of this proposition is based on the following elementary lemma. 

Lemma 3 Let V be a (not necessarily planar) n-gon with external angles 
Pi. Choose a point P and connect it to all vertices ofV. Let ai he the angles 
(see Fig\^ of the triangles at the tip P of the obtained pyramid. Then 

n n 

i=l i=l 

and the equality holds if and only if V is planar and convex^ and the vertex 
P lies inside V. 




Figure 2: Proof of Lemma |31 

Proof. Let us denote by 7^ and 6i the angles of the cite-triangles at the 
vertices of V (see FigEI • The claim of Lemma |31 follows from adding over 
all i = 1, . . . , n two obvious (in)equalities 

ft+i > vr - (7j+i + 5i) 
vr-(7i + (5i) = Oj. 

All inequalities become equalities only in the case when V is planar, convex 
and contains P. Note that some of the external angles (3i may vanish. As a 
corollary we obtain a polygonal version of Fenchel's theorem [2j . 

Corollary 4 

n 



^The obtained pyramid is convex in this case. Note that we distinguish convex and 
strictly convex polygons (and pyramides) . Some of the angles (3i of a convex polygon may 
vanish. The corresponding cite-triangles of the pyramid lie in one plane. 



Proof. For a given V choose the point P varying on a straight hne encircled 
by v. There always exist points P so that the star at P is not strictly 
convex, and thus ^ a^ > 27r. 

Proof of Proposition [^ The claim of Proposition EJ is invariant with 
respect to Mobius transformations. Applying a Mobius transformation M 
which maps the vertex v to infinity, M{v) = oo, one observes that all the 
circles passing through v become straight lines and we arrive at the geometry 
shown in Fig|21with P = Af(oo). Now the claim follows immediately from 
Corollary [l] 

Theorem 5 Let S be a simplicial surface. Then 

W{S) > 0, 

and the equality holds if and only if S is a (part of a) convex polyhedron 
inscribed in a sphere. 

Proof. Only the second statement needs to be proven. Due to Proposition 
mthe equality W{S) = implies that all vertices and edges of 5 are convex 
(but not necessarily strictly convex). Deleting the edges which separate 
triangles lying in one plane one obtains a polyhedral surface Sp with circular 
faces and all strictly convex vertices and edges. Proposition [21 implies that 
for every vertex v there exists a sphere S^ with all vertices of the star S{v) 
lying on it. For any edge (f i, ^2) of Sp two neighboring spheres S^^ and S^^ 
share two different circles of their common faces. This implies S^i = 8^2 
and finally the coincidence of all the spheres S^. 

Discrete conformal energy W defined above is a discrete analogue of the 
Willmore energy for smooth surfaces |14j 

>V(5) = 1 f {h- k2fdA = I H'^dA - I KdA. 
4 Js Js Js 

Here dA is the area element, fci, ^2 the principal curvatures, H = -^{ki + ^2) 
the mean curvature, K = k\k2 the Gaussian curvature of the surface. Here 
we prefer a definition for W with a Mobius invariant integrand. It differs 
from the one in the introduction by a topological invariant. 

Let us mention two important properties of the Willmore energy: 

• yV(5) > and yV(5) = if and only if S is the round sphere. 

• yV(5) (and the integrand {ki — k2f'dA) is Mobius invariant [Hll4j 

Whereas the first claim almost immediately follows from the definition, the 
second one is a non-trivial property. We have shown that the same properties 
hold for the discrete energy W; in the discrete case the Mobius invariance is 
built into the definition and the non-negativity of the energy is non-trivial. 



In the same way one can define conformal (Willmore) energy for simpli- 
cial surfaces in Euclidean spaces of liigher dimension and space forms. 

The discrete conformal energy is well defined for polyhedral surfaces with 
circular faces (not necessarily simplicial). 

3 Computation of the Energy 

Consider two triangles with a common edge. Let a, b,c,d G M.^ be their 
other edges oriented as in Fig|21 Identifying vectors in M^ with imaginary 
quaternions Im H one obtaines for the quaternionic product 

ab = — < a,b > +a X b, (1) 

where < a,b > and a x b are the scalar and vector products in M.^. 




Figure 3: Formula for the angle between circumcircles 

Proposition 6 The external angle (3 € [0, vr] between the circumcircles of 
the triangles in Fig\^ is given by one of the equivalent formulas: 

,^, Re g Re abed 

cos(/3) = -- — - = - = 

\ q\ I abed \ 

< a,c >< b,d > — < a,b >< c,d> — <b,c >< d,a > 

|a||6||c||(i| 

Here q = ab^^cd^^ is the cross-ratio of the quadrilateral. 

Proof. Since Re <?, | <? | and (3 are Mobius invariant it is enough to prove 
the first formula for the planar case a, b,c,d € C, mapping all four vertices 
to a plane by a Mobius transformation. In this case q becomes the classical 
complex cross-ratio. Considering the arguments a,b,c,d G C one easily 
arrives at /? = vr — arg q. The second representation follows from the identity 
b~^ = —b/ \ b \ for imaginary quaternions. Finally applying (^ we obtain 

Re abed =< a, b >< c,d> — <axb,cxd >= 
< a,b >< c,d> + <b,c >< d,a > — < a,c >< b,d > . 



4 Minimizing Discrete Conformal Energy 



Similarly to the smooth Willmore functional W, minimizing the discrete 
conformal energy W makes the surface as round as possible. 

Let us denote by S the combinatorial data of S. The simplicial surface 
S is called a geometric realization of the abstract simplicial surface S. 

Definition 7 Critical points of W{S) are called simplicial Willmore sur- 
faces. The conformal (Willmore) energy of an abstract simplicial surface is 
the infimum over all geometric realizations 



W{S) 



inf Wis). 
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Figure 4: Discrete Willmore spheres of "inscribable" {W = 0) and 
inscribable" iyV > 0) type and discrete Boy surface 
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Kevin Bauer implemented the proposed conformal functional with the 
Brakke's evolver [HI and did some minimization numerical experiments. Ex- 
amples of those are presented in Fig0] In the first and the second lines 
the initial configurations and the corresponding Willmore surfaces minimiz- 
ing the conformal energy respectively are shown. Let us call the gradient 
flow of the energy W the discrete Willmore flow. By this flow the energy 
of the first simplicial sphere decreases to zero and the surface evolves into 
a convex polyhedron with all the vertices lying on a sphere. The abstract 
simplicial surface of the central example is different and we obtain a simpli- 
cial Willmore sphere with positive conformal energy. The third example is 



a simplicial projective plane. The initial configuration is made from squares 
divided into triangles (cf. m|). We see that the minimum is close to the 
smooth Boy surface known to minimize ^^ the Willmore energy for projec- 
tive planes. 

The minimization of the conformal energy for simplicial spheres is related 
to a classical result of Steinitz ^21 who has shown that there exist abstract 
simplicial 3-polytopes without geometric realizations all vertices of which 
belong to a sphere. We call these combinatorial types non-inscribable. 

The non-inscribable examples of Steinitz are constructed as follows jS]. 
Let S be an abstract simplicial sphere with vertices colored in black and 
white. Denote the sets of white and black vertices by V^ and Vb respectively, 
V = VwUVb- Assume that the number of black vertices does not exceed the 
number of white vertices, | V^ \>\ Vb |, and there are no edges connecting 
two white vertices and there are edges connecting black vertices. It is easy 
to see that S with these properties cannot be inscribed in a sphere. Indeed, 
assume that we have constructed such inscribed convex polyhedron. Then 
the equality of the intersection angles at both ends of an edge (see left Fig0 
implies 

27r|y;,|> J^/3(e)>27r|K,|, 

eeE 

and the equalities hold only if all edges connect vertices of different color. 
The obtained contradiction to the assumed inequality implies the claim. 

To construct abstract polyhedra with | Ku |^| ^ | and with some edges 
connecting black points, take a polyhedron P whose number of vertices 
does not exceed the number of faces | F \>\ V |. Color all the vertices 
in black, add white vertices at the faces and connect them to all black 
vertices of a face. We obtain a polyhedron with black (original) edges and 
|Ku|=|-f^|^|^|=|^|- The example with minimal possible number of 
vertices | 1^ |= 11 is shown in FigO The starting polyhedron P here are 
two tetrahedra identified along a common face: F = Q,V = b. 

Hodgson, Rivin and Smith [S] found a characterization of inscribable 
combinatorial types, based on a transfer to the Klein model of hyperbolic 
3-space. It is not clear whether there exist non-inscribable examples of non- 
Steinitz type. 

Numerical experiments lead us to the following 

Conjecture 8 The conformal energy of simplicial Willmore spheres is quan- 
tized 

W = 2TrN, Men. 

Note that this claim belongs to differential geometry of discrete surfaces. 
It would be interesting to find a (combinatorial) meaning of the integer A^. 
Compare also with the famous classification of smooth Willmore spheres by 
Bryant |S1 , who has shown that the energy of Willmore spheres is quantized 

W = 47riV, TV G N. 




Figure 5: A discrete Willmore sphere of "non-inscribable" type with 11 
vertices and W = 27r 

The discrete Willmore energy is defined for the ambient spaces (M" or 
S") of any dimension. This leads to combinatorial Willmore energies 

WniS) = ugWiS), SCS\ 

where the infimum is taken over all realizations in the n— dimensional sphere. 
Obviously these numbers build a non- increasing sequence Wn{S) > Wn+i{S) 
which becomes constant for sufficiently large n. 

Complete understanding of non-inscribable simplicial spheres is an inter- 
esting mathematical problem. However the phenomenon of existence of such 
spheres might be seen as a problem in using of the conformal functional for 
applications in computer graphics, such as fairing of surfaces. Fortunately 
the problem disappears just after one refinement step: all simplicial spheres 
become inscribable. Let S be an abstract simplicial sphere. Define its re- 
finement Sr, as follows: split every edge of S in two by putting additional 
vertices and connect these new vertices sharing a face of S by additional 
edges. 

Proposition 9 The refined simplicial sphere Sr is inscribable, and thus 
WiSR) = 0. 

Proof. Koebe's theorem (see, for example, ^^, |21) claims that every 
abstract simplicial sphere S can be realized as a convex polyhedron S all 
edges of which touch a common sphere 5^. Starting with this realization 
S it is easy to construct a geometric realization Sr of the refinement Sr, 
inscribed in 5^. Indeed choose the touching points of the edges of S with S*^ 
as additional vertices of Sr and project the original vertices of S (which lie 
outside of the sphere 5^) to 5^. One obtaines a convex simplicial polyhedron 
Sr inscribed in S'^. 



Another interesting variational problem with the conformal energy is the 
optimization of triangulations of a given simplicial surface. Here one fixes the 
vertices and chooses an equivalent triangulation (abstract simplicial surface 
S) minimizing the conformal functional. The minimum 

W{V) =mmW{S) 

yields "an optimal" triangulation for a given vertex data. In the case of S"^ 
this optimal triangulation is classical. 

Proposition 10 Let S be a simplicial surface with all vertices V on a two 
dimensional sphere S'^. Then W{S) = if and only if it is the Delaunay 
triangulation on the sphere, i.e. S is the boundary of the convex hull ofV. 

In differential geometric applications like numerical minimizing the Will- 
more energy of smooth surfaces (cf. |lUp it is not natural to preserve the 
triangulation by minimizing the energy, and one should also change the 
combinatorial type decreasing the energy. 

Discrete conformal energy W is not only a discrete analogue of the Will- 
more energy. One can show that it approximates the smooth Willmore 
energy although the smooth limit is very sensitive to the refinement method 
and should be chosen in a very special way. A computation which will be 
published elsewhere shows that if one chooses the vertices of a curvature line 
net of a smooth surface S for the vertices of S and triangularizes it then 
W{S) converges to yV(5) by natural refinement. On the other hand the 
infinitesimal equilateral triangular lattice gives in the limit 3/2 times larger 
energy. Possibly the minimization of the discrete Willmore energy with the 
vertices on the smooth surface could be used for computation of the curva- 
ture line net. We are going to investigate this interesting and complicated 
phenomenon. 

5 Bending of Simplicial Surfaces 

An accurate model for bending of discrete surfaces is important for modelling 
in virtual reality. 

Let Sq be a thin shell and S its deformation. The bending energy of 
smooth thin shells is given by the integral 

E= l(H-HnfdA, 



where Hq and H are the mean curvatures of the original and deformed 
surface respectively. For Hq = it reduces to the Willmore energy. 
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To derive the bending energy for simplicial surfaces let us consider the 
limit of fine triangulation, i.e. of small angles between the normals of neigh- 
boring triangles. Consider an isometric deformation of two adjacent tri- 
angles. Let 9 be the complement of the dihedral angle of the edge e, or, 
equivalently, the angle between the normals of these triangles (see FigEI 
and P{9) the external intersection angle between the circumcircles of the 
triangles (see Fig^ as a function of 9. 

Proposition 11 Assume that the circumcenters of the circumcircles of two 
adjacent triangles do not coincide. Then in the limit of small angles 9^0 
the angle jS between the circles behaves as follows: 

(5{9)=(5{<d)+^-9^ + o{9^). 
ij 

Here I is the length of the edge and L ^ is the distance between the centers 

of the circles. 

This proposition and our definition of conformal energy for simplicial 
surfaces motivate to suggest 

I 



^ = E 



l'' 



ee_B 

for the bending energy of discrete thin-shells. 

/ 




Figure 6: To definition of the bending energy for simplicial surfaces 

In ^ [7j similar representations for the bending energy of simplicial sur- 
faces were found empirically. They were demonstrated to give convincing 
simulations and good comparison with real processes. In [7j the distance 
between the barycenters is used for L in the energy expression but possible 
numerical advantages in using circumcenters are indicated. 

Using the Willmore energy and Willmore flow is a hot topic in computer 
graphics. Applications include fairing of surfaces and surface restoration. 
We hope that our conformal energy will be useful for these applications and 
plan to work on them. 
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